Award  Number:  W81XWH-04-1-0239 


TITLE:  Early  Detection  of  Breast  Cancer  by  Fluorescence  Molecular  Tomography 


PRINCIPAL  INVESTIGATOR:  Vasilis  Ntziachristos,  Ph.D. 

Stephen  D.  Windsor 


CONTRACTING  ORGANIZATION:  The  General  Hospital  Corporation 

Massachusetts  General  Hospital 
Boston,  MA  02114 


REPORT  DATE:  July  2005 


TYPE  OF  REPORT:  Annual 


PREPARED  FOR:  U.S.  Army  Medical  Research  and  Materiel  Command 
Fort  Detrick,  Maryland  21702-5012 


DISTRIBUTION  STATEMENT:  Approved  for  Public  Release; 

Distribution  Unlimited 


The  views,  opinions  and/or  findings  contained  in  this  report  are  those  of  the  author(s)  and 
should  not  be  construed  as  an  official  Department  of  the  Army  position,  policy  or  decision 
unless  so  designated  by  other  documentation. 


20060503051 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and  maintaining  the 
data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for  reducing 
this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202- 
4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently 
valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


2.  REPORT  TYPE  3.  DATES  COVERED  (From  -  To) 

Annual  1  Jul  04  -  30  Jun  05 


1.  REPORT  DATE  (DD-MM-YYYY) 

01-07-2005 


4.  TITLE  AND  SUBTITLE 

Early  Detection  of  Breast  Cancer  by  Fluorescence  Molecular  Tomography 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 

W81XWH-04-1-0239 


5c.  PROGRAM  ELEMENT  NUMBER 


6.  AUTHOR(S) 

Vasilis  Ntziachristos,  Ph.D. 
Stephen  D.  Windsor 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


E-Mail:  vasilis@helix.mgh.harvard.edu 


5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

The  General  Hospital  Corporation 
Massachusetts  General  Hospital 
Boston,  MA  02114 


8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 


9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

U.S.  Army  Medical  Research  and  Materiel  Command 
Fort  Detrick,  Maryland  21702-5012 


10.  SPONSOR/MONITOR’S  ACRONYM(S) 


11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  Public  Release;  Distribution  Unlimited 


14.  ABSTRACT 

Early  detection  of  breast  tumors  continues  to  be  a  significant  diagnostic  challenge  and  thus 
remains  the  focus  of  attention  of  a  number  of  medical  research  groups.  Molecular  targeting 
approaches  have  tremendous  potential  for  early-detection  because  they  rely  on  elucidation  of 
abnormal  gene-expression,  rather  than  on  discovery  of  retarded  anatomical  changes  inflicted 
by  growing  tumors  upon  their  microenvironment.  Originally,  we  proposed  an  investigation  of 
whether  highly-sensitive  fluorescence  molecular  tomography  (FMT)  could  be  used  to  detect 
breast  cancer  at  its  earliest  and  smallest  stages  via  the  detection  of  injected,  protease- 
activatable  molecular  probes.  We  enter  year  2  with  the  year  1  goals  fully  achieved,  having 
completed  construction  of  the  proposed  FMT  imaging  system  and  having  optimized  the  system 
using  phantoms.  Furthermore,  we  have  made  significant  progress  toward  accomplishing  the  goals 
for  year  2  by  beginning  the  in-vivo  early  detection  and  screening  studies  ahead  of  schedule. 


15.  SUBJECT  TERMS 

Molecular  imaging,  fluorescence,  tomography,  early  detection,  specificity,  breast  cancer 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

OF  ABSTRACT 

18.  NUMBER 

OF  PAGES 

19a.  NAME  OF  RESPONSIBLE  PERSON 

USAMRMC 

a.  REPORT 

U 

b.  ABSTRACT 

U 

c.  THIS  PAGE 

u 

UU 

31 

19b.  TELEPHONE  NUMBER  (include  area 
code) 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  Z39.18 


Table  of  Contents 


Cover . 1 

SF  298 . 2 

Introduction . 4 

Body . 7 

Key  Research  Accomplishments . 18 

Reportable  Outcomes . 18 

Conclusions . 19 

References . 19 


Appendices 


22 


< 


Introduction 

Breast  cancer  is  the  second  most  common  female  malignancy  in  the  United  States.  The  ability  to 
detect  malignancies  in  their  earliest  stages  is  of  tremendous  importance.  Early  detection  has  been 
shown  to  increase  likelihood  of  survival  [1,  2]  and  to  reduce  the  risk  of  disease  recurrence  [3  - 
6].  However,  early  detection  continues  to  be  a  significant  diagnostic  challenge  and  thus  remains 
the  focus  of  attention  of  a  number  of  medical  research  groups. 

Well-established  detection  techniques,  most  notably  mammography,  remain  limited  in  their 
sensitivity  to  small-volume  disease.  Though  mammography  has  been  of  great  importance  for 
detection  of  breast  malignancies,  many  tumors  have  evaded  detection  by  this  method  largely  due 
to  its  limited  contrast.  The  technique  is  particularly  hampered  by  background  interference  such 
as  increased  tissue  density  or  scarring  from  prior  surgery.  Functional  imaging,  which  includes 
Magnetic  Resonance  Imaging  (MRI)  and  Positron  Emission  Tomography  (PET),  has  been  used 
as  a  compliment  to  anatomical  imaging  and  has  shown  significant  potential  to  image  functional 
characteristics  of  tumors  [7],  However,  the  high  cost  and  complex  operational  logistics  of  these 
modalities  may  limit  their  practicality  for  use  in  large-scale  screening  studies. 

An  additional  parameter  that  plays  a  vital  role  in  breast  cancer  detection  is  the  ability  to  detect 
specific  disease  signatures  based  on  the  molecular  profile  of  the  disease.  Molecular  imaging,  i.e., 
imaging  of  cellular'  and  sub-cellular  processes  that  are  associated  with  molecular  onset  and 
evolution  of  disease,  can  lead  to  increased  sensitivity  and  specificity  and  can  improve  the  overall 
performance  of  a  diagnostic  study.  Improved  specificity,  based  on  the  detection  of  molecular 
profiles,  can  lead  to  the  overall  reduction  in  unnecessary  biopsies,  which  are  associated  with 
increased  healthcare  costs  and  patient  distress,  and  thus  can  lead  to  personalized  medicine.  It  is 
foreseen  that  such  methods  will  help  optimize  therapeutic  strategies  and  will  enable  the  effects  of 
treatment  to  be  more  closely  followed  in  time  scales  of  hours  or  days  instead  of  the  longer 
timeframe  of  weeks  or  months  that  is  the  current  standard  for  on  anatomical  imaging  methods. 
Some  of  our  recent  review  articles  [8,  9]  provide  more  detailed  discussions  of  molecular  imaging 
in  general,  and  one  of  our  recent  studies  [10]  demonstrates  the  ability  of  optical  imaging  to 
monitor  disease  behavior. 

Optical  methods  have  emerged  as  economic,  safe,  and  effective  alternatives  for  functional 
imaging.  Diffuse  Optical  Tomography  (DOT)  has  been  developed  for  imaging  of  breast  cancers 
because  of  its  ability  to  non-invasively  quantify  oxy-  and  deoxy-hemoglobin  concentrations, 
which  can  be  used  to  characterized  angiogenesis  and  hypoxia  [1 1-15].  In  addition,  DOT  has  been 
used  in  combination  with  MRI  to  detect  tumors  based  on  the  extrinsically-administered  contrast 
agent  Indocyanine  Green  (ICG)  that  marks  angiogenesis  and  permeability.  However,  these 
markers,  though  they  enable  DOT  to  effectively  image  well-developed  tumors,  are  not  ideal  for 
detection  of  the  earliest  stages  of  disease.  Consequently,  while  DOT  may  certainly  be  a  useful 
tool  for  evaluating  larger  tumors,  a  more  effective  optical  modality  is  needed  for  early  detection. 
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Fluorescence  Molecular  Tomography  (FMT)  is  a  relatively  new  optical  tomographic  imaging 
approach  that  capitalizes  on  the  strengths  of  molecular  imaging.  It  bases  detection  on  specific 
molecular  signatures,  using  highly  targeted  fluorescent  probes,  rather  than  on  anatomical  or 
functional  changes  within  tumor  microenvironments.  This  approach  yields  a  shift  in  the 
radiological  paradigm  that  has  traditionally  focused  on  anatomy  and  physiology  [16J.  Unlike 
DOT,  FMT  can  three-dimensionally  resolve  protein  expression,  protease  activity,  receptor 
regulation,  and  similar  molecular  markers  that  play  vital  roles  in  carcinogenesis.  It  has  been 
demonstrated  in-vivo  that  NIR  photons  can  penetrate  more  than  12  centimeters  into  the  breast, 
and  that  fluorochromes  accumulating  in  small  tumors  are  detectable  through  the  human  breast 
[17].  We  have  previously  shown  that  FMT  is  capable  of  resolving  phantoms  placed  less  than  1 
mm  apart  [18],  and  the  sensitivity  of  this  technology  has  been  predicted  to  be  adequate  for 
detection  of  human  breast  lesions  smaller  than  5  mm  [19].  More  recently  it  has  been  shown  that 
the  technique  can  visualize  apoptosis  as  a  response  to  treatment  [10].  Furthermore,  because  non¬ 
ionizing  radiation  is  used,  the  technology  allows  for  repeated  imaging  for  screening  or 
monitoring  of  disease.  FMT  thus  holds  great  promise  for  early  detection  of  breast  malignancies. 

At  the  time  of  our  grant  proposal,  several  elegant  targeted  fluorescent  probes  had  already  been 
developed  [20-22]  for  use  with  optical  methods  such  as  FMT.  Among  the  most  notable  of  these 
were  the  so-called  “smart”  probes  (also  termed  “activatable  probes”),  which  are  optically  silent 
until  they  interact  with  a  specific  enzyme,  such  as  a  protease,  which  then  activates  the  probe, 
converting  it  to  a  brightly  fluorescent  state  [20].  These  probes  are  powerful  not  only  because  of 
their  specificity  for  cancer- associated  enzymatic  activity,  but  also  because  they  significantly 
minimize  background  levels  and  exhibit  several-fold  amplification  of  fluorescence  upon 
activation,  enabling  dramatic  improvements  in  signal-to-noise  contrast.  Because  activatable 
probes  fluorescence  only  in  the  presence  of  specific  tumor  phenotypes,  they  can  easily  be 
localized  and  quantified  against  a  virtually  “dark”  background,  allowing  for  highly  sensitive  and 
specific  detection. 

Our  proposal  described  using  a  smart  probe  that  is  activated  by  major  cathepsins,  which  are 
members  of  the  cysteine  family  of  proteases  and  play  a  central  role  in  tumorigenesis  and  matrix 
invasion  [23-26].  For  example,  it  has  been  demonstrated  that  high  levels  of  cathepsin  B 
expression  correlates  positively  with  aggressive  behavior  and  progression  of  human  tumors  and 
negatively  with  patient  survival  [23],  We  therefore  chose  to  evaluate  the  performance  of  this 
probe  for  imaging  breast  cancer  development  and  progression. 

Until  the  study  was  supported  by  this  grant,  in-vivo  imaging  of  such  probes  were  limited  to  basic 
feasibility  studies,  generally  involving  imaging  of  xenographic  implants  in  nude  mice.  Generally, 
it  was  unknown  how  early  in  disease  progression  such  probes  could  be  used  to  reliably  detect 
malignancies,  as  the  in-vivo  sensitivity  limits  of  FMT  have  not  been  rigorously  explored.  The 
goals  of  our  proposal  were  primarily  (1)  to  develop  and  optimize  optical  technology  for  highly 


5 


sensitive  detection;  and  (2)  to  evaluate  the  early- detection  limits  of  FMT  by  conducting  a  blind 
screening  study  of  transgenic  mice  that  develop  spontaneous  tumors  of  the  mammary  fat  pad, 
closely  mimicking  the  natural  development  and  progression  of  human  breast  cancer.  We 
proposed  developing  imaging  tools  designed  to  systematically  reduce  background  noise  levels 
and  to  improve  fluorochrome  localization  and  quantification.  The  study  was  designed  to 
investigate  screening  efficacy  and  to  predict  sensitivity  limits  of  human  breast  cancer. 

The  hypothesis  of  this  proposal  is  that  highly  sensitive  fluorescence  imaging  (tomosynthesis)  of 
molecular  contrast  will  enable  earlier  detection  of  breast  tumors  based  on  signals  associated  with 
tumor  growth  and  matrix  invasion.  We  predicted  that  this  technology  could  achieve  high 
detection  specificity. 

Our  study  had  three  major  goals: 

Aim  1.  To  develop  appropriate  photon  illumination  and  detection  strategies  and  reconstruction 
methods  for  highly  sensitive  fluorescence  detection  and  quantification. 

Aim  2.  To  utilize  optimal  photon  technology  from  Aim  1  to  study  in-vivo  detection  of 
spontaneous  disease  in  transgenic  mice  and  to  follow  local  disease  progression  and  distant 
metastases. 

Aim  3.  To  translate  the  animal  study  findings  into  clinically-relevant  detection  schemes. 

The  first  of  these  aims  was  designed  to  find  an  optimal  illumination  technology  and  processing 
method  in  order  to  obtain  the  highest  detection  sensitivity  of  fluorochromes  in  tissues.  We 
proposed  investigating  limited  angle  projections  versus  multi-angle  illumination,  planar 
illumination,  and  transillumination  schemes. 

The  second  aim  focused  on  the  clinical  relevance  of  FMT.  We  proposed  a  screening  study 
involving  transgenic  mice  that  spontaneously  develop  breast  tumors.  Our  interest  was  and 
continues  to  be  primarily  (1)  to  characterize  detection  capacity  as  a  function  of  disease 
progression  and  tumor  growth;  and  (2)  to  quantify  fluorescence  concentration  in-vivo  and 
determine  whether  it  correlates  with  tumor  expression  profiles.  This  screening  study  would 
enable  us  to  establish  early-detection  limits. 

The  third  aim  was  to  translate  the  copious  data  acquired  from  exploring  the  first  two  aims  into  an 
evaluation  of  the  clinical  feasibility  of  FMT  for  early-detection  of  breast  cancer. 

In  the  grant  proposal  we  presented  a  timeline  for  accomplishing  our  primary  aims  for  years  1-3 
of  the  study  (Table  1). 
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Table.  1  Overview  of  aims  and  time  line 


Year  01 

Year  02 

Year  03 

Aim  1:  Method  Development  &  Optimization 

•  Construction  of  4-mode  imaging  chamber 

x 

•  Developing  imaging  algorithms 

X 

•  Optimizing  imaging  parameters  with  phantoms 

•  Optimize  imaging  parameters  with  orthotopic  mouse  model(s) 

X 

X 

Aim  2:  In  vivo  Screeninq  study 

•  Establish  appropriate  transgenic  mouse  model(s) 

•  Breed  and  screen  mice  for  early  detection  of  tumors 

x 

X 

X 

■  Study  local  and  distant  disease  progression. 

x 

X 

Aim  3:  Clinical  translation 

•  Construction  of  breast-like  phantoms 

x 

X 

•  Experiments  with  phantoms 

X 

♦  Simulations 

X 

To  date,  we  have  not  only  met,  but  also  exceeded  the  goals  for  year  1,  and  we  have  also  made 
significant  progress  in  reaching  the  goals  for  year  2.  The  following  sections  provided  a  detailed 
evaluation  of  our  progress. 

Body 

We  have  completed  all  of  the  goals  set  for  year  1,  which  encompasses  all  of  Aim  1  and  a  small 
portion  of  Aim  2.  Furthermore,  we  have  made  significant  progress  towards  achieving  Aim  2  and 
thus  have  a  strong  head-start  in  completing  the  goals  for  year  2.  The  following  is  a  detailed 
description  of  accomplishments  and  reported  outcomes  achieved  over  the  past  year: 

Aim  1.  To  develop  appropriate  photon  illumination  and  detection  strategies  and 
reconstruction  methods  for  highly  sensitive  fluorescence  detection  and  quantification. 


Al.l:  Construction  of  an  integrated  illumination/detection  small  animal  imaging  system. 

We  constructed  a  4-mode  imaging  chamber,  a  schematic  of  which  is  presented  in  Figure  1. 
Photographs  of  the  system  are  shown  in  Figure  2.  This  system  is  based  on  a  previously-reported 
scanner  [18].  Illumination  is  provided  by  either  a  672  nm  wavelength  or  a  748  nm  wavelength 
cw  laser  (a)  routed  to  a  two-channel  optical  switch  (b).  One  output  of  this  switch  connects  to  an 
illumination  branch  (f)  for  reflectance  imaging  (g),  while  the  other  output  connects  a  source  fiber 
(c)  to  a  1x46  optical  switch  (d)  which  establishes  all  the  other  three  imaging  modes.  The  1x46 
switch  multiplexes  the  single  light  input  (c)  to  46  output  fibers  (e),  which  couple  the  laser  light 
from  the  switch  to  the  animal  chamber  (h).  The  46  fibers  are  inserted  through  the  backside  of  the 
animal  chamber,  abutting  the  inner  back  wall,  and  are  arranged  equidistantly  in  a  hexagonal 
mesh.  The  chamber  is  identical  to  ones  that  have  been  described  previously  (Graves  et  al  2003). 
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Furthermore,  since  the  grant  proposal  we  have  constructed  a  rotational  stage  and  mouse  holder 
(Figure  2,  right)  that  may  be  inserted  into  the  chamber  to  achieve  complete  projection  (360) 
geometry.  A  low  noise,  512x512  element  cooled  CCD  camera  (k)  collects  photons  that  pass 
through  appropriate  bandpass  filters  (i).  Data  is  then  processed  using  a  3.20  GHz  Pentium  4 
personal  computer  with  2  GB  of  RAM  running  Windows  XP  (n).  The  system  accommodates  all 
four  illumination  schemes  that  were  described  in  the  grant  proposal,  namely  reflectance  imaging, 
transillumination,  tomosynthesis,  and  tomography.  These  illumination  schemes  are  also  shown  in 
Figure  3. 


Figure  1:  FMT  imaging  system. 


(middle),  and  a  photograph  of  the  rotational  stage  with  a  mouse  in  place  in  the  holder  (right), 
which  can  be  placed  into  the  imaging  chamber  for  complete  projection  (360°)  geometry. 
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(A)  Planar  (reflectance)  (B,  iranslllumlnatlon 


Figure  3:.  Schematics  of  the  4  modes  of  imaging 


Reflectance  Imaging 

The  illumination  branch  (Figure  1,  f,  g;  Figure  3,  middle)  illuminates  the  animal  through  the 
chamber  window,  and  reflected  intrinsic  or  fluorescence  photons  are  recorded  on  the  CCD 


Figure  4:  Intrinsic  (left)  and  fluorescence  (right)  reflectance  images  of  a  mouse  with 


subcutaneous  tumors  (red  arrows  in  right  image).  The  liver  is  indicated  with  a  white 


arrow.  The  mouse  was  injected  with  2  nmol  of  a  cathepsin-activatable  probe  24  hours 


prior  to  imaging. 
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camera  using  band-pass  filters.  Figure  4  presents  intrinsic  and  fluorescence  reflectance  images  of 
a  mouse  with  two  subdermal  tumors  within  the  mammary  fat  pads,  the  result  of  subcutaneous 
implantation  of  3  x  106  HT-1080  cells  into  each  tumor  site.  Twenty- four  hours  prior  to  imaging, 
the  mouse  was  injected  with  2  nmol  of  a  cathepsin-activatable  NIR  probe.  The  intrinsic  image  is 
useful  for  anatomical  coregistration  of  the  fluorescence  image,  which  clearly  shows  two  tumors 
as  well  as  the  liver,  which  is  known  to  activate  the  probe. 

T ransillumination . 

Each  of  the  46  source  fibers  directed  toward  the  back  of  the  imaging  chamber  (Figure  1,  e)  is 
sequentially  switched  on,  the  transmitted  pattern  is  filtered  with  appropriate  intrinsic  or 
fluorescence  band  pass  filters,  and  photons  are  recorded  on  the  CCD  camera.  The  resulting 
images  are  added  together  to  create  a  composite  image,  which  is  geometrically  similar  to  an  X- 
ray  mammogram.  The  fluorescence  composite  image  may  also  be  divided  by  the  intrinsic 
composite  image  to  generate  a  normalized  image,  which  significantly  reduces  sensitivity  to 
background  heterogeneity  and  thus  provides  a  more  accurate  image.  Another  effective  approach 
is  to  normalize  each  of  the  46  fluorescence  images  individually  by  their  corresponding  intrinsic 
images  prior  to  summation  (Figure  5).  Figure  6  demonstrates  an  example  of  transillumination 
imaging  on  a  mouse  with  a  mammary  tumor  implanted  in  the  right  side  of  the  animal.  The 


Figure  5:  Construction  of  a  transillumination  composite  image  of  a  mouse 
with  subcutaneous  mammary  tumors;  the  images  are  of  a  region  of  interest 
chosen  to  encapsulate  the  tumors.  The  mouse  was  injected  with  a  cathepsin- 
activatable  NIR  probe  24-hours  prior  to  imaging.  Fluorescence 
transillumination  images  (A)  were  divided  by  corresponding  intrinsic 
transillumination  images  (B).  Low  thresholds  were  applied  to  the  images 
prior  to  normalization  to  reduce  background  noise.  The  individually- 
normalized  images  were  added  to  create  a  composite  image  (C)  that 
provides  a  more  accurate  localization  of  fluorescence  signal.  Arrowheads 
indicate  tumor  locations. 
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Fluorescence 

Photograph  Transillumination 


Figure  6:  Transillumination  imaging  of  breast  tumors,  (a)  Photograph  of  the  animal. 

The  arrow  indicates  the  position  of  the  tumor,  (b)  Inverted  transillumination  image 
after  data  normalization  reveals  higher  fluorescence  from  the  tumor  (darker  indicates 
higher  signal),  which  is  associated  with  up-regulation  of  cathepsins. 

animal  was  injected  with  a  protease  activated  probe  24-hours  prior  to  imaging.  Generally, 
transillumination  imaging  is  better  suited  for  volumetric  measurements  compared  to  reflectance 
imaging  because  the  whole  animal  volume  is  probed.  This  is  because  light  is  passing  through  the 
animal  and,  therefore,  excites  all  fluorescence  present  in  its  path. 

Tomosynthesis. 

We  had  made  significant  progress  with  implementing  limited  angle  projection  FMT,  the 
equivalent  of  FMT  tomosynthesis.  For  this  approach  data  is  gathered  in  the  same  manner  as  in 
the  transillumination  approach;  however,  unlike  that  approach,  tomosynthesis  involves  applying 
algorithms  to  the  data  that  solve  a  volume  integral  equation  relating  spatially-distributed 
fluorochrome  concentrations  (in  physical  units)  to  the  transillumination  measurements.  Inversion 
is  performed  using  an  algebraic  reconstruction  technique  with  convergence  assumed  after  20 
iterations,  as  empirically-determined  to  yield  optimal  imaging  performed  based  on  phantom 
measurements.  This  process  yields  a  three-dimensional  “map”  of  fluorochromes  concentration 
that  is  subsequently  sliced  for  viewing  (Figure  7). 


z=T,aacm  z=1.1<Scm 


Figure  7:  Tomographic  slices  of  a  mouse  with  subcutaneous  tumors. 
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Tomography. 

To  achieve  superior  imaging  performance  it  is  important  to  illuminate  using  a  large  number  of 
projections  and  detect  signal  around  the  animal,  similarly  to  other  tomographic  techniques  such 
as  X-ray  CT,  PET  or  SPECT.  For  this  reason  we  have  implemented  the  rotational  device  (Figure 
1,  m)  which  rotates  the  object  of  interest  in  front  of  the  illumination  path  and  CCD  camera  to 
implement  a  number  of  projections.  The  system  attains  rotational  accuracy  of  better  than  0.5 
degree,  but  generally  a  fewer  number  of  projections  is  used  (typically  varying  from  72  to  18). 
Figure  8  shows  an  example  of  imaging,  in  this  case  also  utilizing  photon  pulses  and  time-gated 
detection.  In  this  example  a  solid  phantom  was  utilized,  made  of  polyester  resin  in  which  Ti02 
particles  and  India  Ink  were  added  to  yield  absorption  and  scattering.  The  phantom  was  cut  in  the 
cross-sectional  shape  of  a  rectangle  and  a  square  as  shown  in  Figure  8a  and  was  immersed  in  a 
1%  intralipid  solution  which  similarly  also  contained  India  Ink  to  yield  absorption.  The  optical 
properties  of  the  intralipid-ink  solution  were  (|la=0.1  cm'1  (4’=  7  cm'1)  whereas  the  solid 
phantom  demonstrated  8x  this  background  attenuation. 


Figure  8:  Complete  projection  tomography  imaging  here  in  a  diffusive  phantom  of  a 
triangle  and  a  square,  (a)  Schematic  of  the  arrangement  of  the  phantom,  (b)  360° 
projection  reconstructed  image  using  early  photons,  (c)  360°  projection  reconstructed 
image  using  CW. 


The  reconstructed  images  are  shown  in  Figures  8b  and  c  respectively.  Figure  8b  demonstrates  the 
ability  of  early  photons  to  reconstruct  correctly  the  position  and  shape  of  the  phantoms.  Figure 
8c,  which  is  based  on  CW  data  also  reconstructs  the  position  of  the  phantoms  but  cannot 
reconstruct  shape  due  to  the  diffusive  nature  of  the  phantom.  However,  although  contrast  is 
reduced,  the  position  and  size  of  the  phantoms  is  also  reconstructed.  In  year  2  we  will  continue 
optimizing  the  parameters  for  reconstruction  of  CW  data  based  on  complete  projection 
tomography  data  and  expect  superior  reconstruction  to  be  achieved. 


A  1.2:  Development  of  imaging  software 

In  this  first  year  we  have  achieved  the  following  advancements  in  terms  of  imaging  software: 

1.  Developed  algorithms  for  transillumination  imaging  using  appropriate  normalized 
techniques  to  improve  contrast  and  robustness  (1  publication  in  press  to  appear  in  J.  of 
Biomedical  Optics) 

2.  Developed  advanced  background  fluorescence  techniques  based  on  physical  models  and 
improve  on  tomographic  performance  and  contrast  (1  publication  in  preparation  ) 

3.  Evaluated  algorithms  for  robust  imaging  in  highly  heterogeneous  backgrounds  as 
appropriate  for  in-vivo  imaging  (One  publication  in  press  to  appear  in  IEEE  transaction 
on  Medical  Imaging) 

4.  Adapted  algorithms  already  developed  for  tomosynthesis  to  complete-projection  imaging. 

A1.3:  Phantom  studies 


We  performed  several  phantom  studies  to  evaluate  the  different  algorithms  developed  and  assess 
the  overall  performance  of  the  method.  In  particular  we  constructed  highly  heterogeneous 
phantoms  using  combinations  of  solid  and  liquid  diffusive  media  and  investigated  the 
performance  of  the  methods  with  a  varying  degree  of  heterogeneity  (for  examples  see  appendix). 
In  addition  we  streamlined  the  developed  of  resin  phantoms  (Figure  9),  of  different  optical 
properties  that  can  simulate  tissues  of  different  optical  properties  (for  example  varying  optical 
properties  in  the  breast  due  to  aging)..  We  have  further  examined  the  effects  of  background 
fluorescence  in  homogeneous  and  heterogeneous  phantoms  (paper  in  preparation).  We  further 
confirmed  previous  studies  [19,  27]  with  liquid  phantoms  attaining  the  optical  properties  of  the 
human  breast  that  detection  sensitivity  is  of  the  order  of  10  nM  (volume:  200  (il)  of  fluorescence 
dye  (Cy5.5  ex  wvl:  672  nm,  max  em  wvl  690  nm;  5  sec  exposure  per  source)  at  a  phantom  with 
the  optical  properties  of  the  breast  at  2  cm  width.  In  accordance  to  our  time-plan,  in  year  2  we 
will  build  breast  like  phantoms  and  test  sensitivity  limits  with  phantoms  that  bear  geometrical 
resemblance  to  the  breast  to  further  evaluate  these  figures. 


Figure  9:  Two  photographs  of  resin  phantoms  of  different 
shapes  and  varying  degrees  of  absorbance  and  fluorescence. 
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A1.4:  Evaluation  of  imaging  practices. 


We  further  undertook  studies  to  compare  the  performance  of  planar  and  tomographic  methods  to 
quantify  fluorochromes  as  a  function  of  depth  and  optical  properties.  Here  we  summarize  some 
of  the  findings,  although  more  details  can  be  found  in  the  appendix. 


The  first  study  compared  planar  imaging  and  tomosynthesis  (tomographic)  performance  in  a  side 
by  side  comparison.  Shown  in  FigurelOa  is  imaging  of  two  ~1.5mm  diameter  fluorescent  tubes 
(500nM  Cy5.5,  3mm  apart)  immersed  at  different  depths  in  a  water  solution  of  intralipid  and 
India  ink  simulating  the  optical  properties  of  tissue.  Both  modalities  resolve  the  tubes  when 
placed  superficially  although  FMT  offers  better  resolution.  However,  detection  becomes  highly 


Figure  10:  Side-by-side  comparison  of  planar  and  tomographic  imaging 
(a)  as  a  function  pf  depth,  (b)  as  a  function  of  varying  local  optical 
properties,  and  (c)  as  a  function  of  heterogeneity.  In  all  cases 
tomosynthesis  (tomography)  is  shown  superior  to  planar  imaging. 


challenging  as  a  function  of  depth  for  planar  imaging  with  the  tubes  barely  visible  even  3mm 
deep.  Depth  sensitivity  highly  depends  on  the  size  and  strength  of  fluorochrome,  but  the  images 
contrast  the  superior  ability  of  tomography  to  look  deeper  in  diffuse  media  and  with  higher 
resolution,  using  in  this  case  identical  hardware  to  the  one  used  for  planar  imaging.  The  color 
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images  are  reconstructed  slices  obtained  at  different  depths  and  superimposed  on  gray  scale 
photographs  of  the  tubes  for  visualization  purposes.  Similar  results  are  obtained  in  comparisons 
of  transillumination  and  tomography,  although  transillumination  can  still  record  objects  deeper 
than  in  planar  imaging.  Similarly,  the  ability  to  image  fluorochromes  with  varying  background 
optical  properties  is  showcased  in  Figure  10b,  which  depicts  planar  and  reconstructed  images  of 
two  fluorescent  tubes,  the  left  tube  containing  double  the  amount  of  Cy5.5  in  the  right  tube 
(400nM  vs.  200nM).  The  tubes  are  immersed  in  the  same  diffuse  fluid  as  in  (a).  Both  techniques 
accurately  resolve  the  2:1  relation  in  fluorochrome  concentration  between  left  and  right  tubes 
when  background  absorption  is  the  same  in  both  tubes  (top  row).  However  when  India  Ink  is 
added  in  the  left  tube  to  simulate  a  three-fold  increase  in  vascularization  (absorption),  planar 
performance  is  significantly  reduced,  erroneously  reporting  a  1:1  Cy5.5  concentration  in  the  two 
tubes  (bottom  row).  This  is  because  the  added  ink  absorbs  more  fluorescent  photons.  In  contrast 
tomography  can  correct  for  these  phenomena  and  demonstrates  more  robust  performance 
reporting  1.8:1  relation  in  this  case.  Finally,  Figure  10c  demonstrates  the  capacity  of  tomography 
to  resolve  a  single  fluorescent  tube  in  the  middle  of  a  highly  diffusive  medium  when  it  is 
asymmetrically  surrounded  by  5  strong  absorbers  of  2x  and  3x  the  background  absorption  and  of 
comparable  sizes.  This  result  demonstrates  the  robustness  of  tomography  to  image  highly 
heterogeneous  media,  an  ability  that  is  mainly  due  to  the  data  normalization  techniques  (image 
labeled  nbom)  employed  tomographically.  Absence  of  normalization  results  in  distorted 
tomographic  imaging  performance  (shown  in  the  image  labeled  nBom).  Both  images  are  slices 
from  a  three  dimensional  reconstruction  obtained  at  9mm  depth.  Planar  imaging  could  not  detect 
the  presence  of  a  fluorochrome  in  this  highly  absorbing  and  heterogeneous  background. 

Overall,  these  findings  exemplify  that  planar  imaging  should  be  used  with  caution.  Signal 
intensity  relates  linearly  to  fluorochrome  concentration  but  non-linearly  to  depth,  size  and  optical 
properties  and  is  further  complicated  by  the  high  scattering  nature  of  tissue.  Tomosynthesis  and 
tomography  have  the  potential  to  improve  on  these  limitations  and  to  offer  more  robust  and 
accurate  imaging. 

Aim  2.  To  utilize  the  optimal  photon  technology  found  in  Aiml  to  study  in  vivo  detection  of 
spontaneous  disease  in  transgenic  mice  and  to  follow  local  disease  progression,  distant 
metastases  and  treatment  response. 

As  stated  in  the  grant  proposal,  we  selected  a  transgenic  mouse  model  (c-neu  oncogene;  strain 
name  FVB/N-TgN(MMTVneu)202Mul,  Jackson  Laboratory,  Bar  Harbor,  Maine)  that  accurately 
mimics  known  genetic  alterations  seen  in  human  breast  cancer.  In  this  model  the  mouse 
mammary  virus  tumor  (MMTV)  promoter  directs  expression  of  the  c-neu  oncogene  to  the 
mammary  fat  pads.  The  model  was  originally  reported  to  develop  mammary  tumors  in  50%  of 
females  by  89  days,  but  ultimately  there  is  90%  penetrance  by  5-6  months.  The  model  is  of  high 
biological  relevance  because  the  Her-2/neu  gene  is  over-expressed  in  over  90%  of  human  breast 
cancers. 
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A2.3:  Imaging  studies 


Thus  far  16  mice  have  resulting  from  breeding  and  have  completed  screening.  Currently  an 
additional  16  mice  are  approaching  the  age  of  tumor  onset  and  will  soon  be  enter  imaging 
studies.  Therefore,  we  foresee  successful  completion  of  imaging  studies  involving  a  total  of  60 
transgenic  mice  (40  experimental  mice,  10  control  mice,  and  10  mice  for  purposes  of  parameter 
optimizations)  within  the  allotted  time  as  originally  stated  in  the  grant  proposal. 

As  stated  previously  in  the  grant  proposal,  the  purpose  of  the  imaging  study  is  to  explore  the 
early-detection  and  disease-monitoring  abilities  of  FMT.  Figure  1 1  presents  results  from  a  mouse 
with  a  developing  tumor  and  showcases  the  ability  of  the  system  to  follow  tumor  development 
over  time.  Furthermore,  it  was  found  that  fluorochrome  concentration,  as  calculated  by  FMT, 
increased  monotonically  with  tumor  diameter  (Figure  12),  an  expected  and  important  result,  as  it 
is  predicted  that  probe  will  accumulate  and  activate  selectively  within  tumors  and  that, 
consequently,  larger  tumors  should  display  higher  concentration  of  activated  probe. 


Figure  11:  A  mouse  imaged  prior  to  visible  manifestation  of  tumors  and  subsequently  imaged  on 
a  weekly  basis  to  follow  tumor  growth  from  3  mm  to  6  mm  in  diameter.  The  mouse  was  injected 
via  tail-vein  with  2  nmol  of  a  cathepsin-activatable  NIR  probe  24-hours  prior  to  each  imaging 
session. 


16 


Fluorochrome  Concentration  (jjM) 


Figure  12:  Fluorochrome  concentration  from  the 
mouse  shown  in  Figure  1 1  as  a  function  of  tumor 
diameter  as  measured  with  calipers. 

Using  a  different  probe,  consisting  of  the  Her-2  antibody  conjugated  to  Cy5.5,  also  injected  into 
a  transgenic  mouse  24-hours  prior  to  imaging,  the  system  was  able  to  resolve  a  1  mm  tumor  (as 
measured  across  the  diameter  of  the  tumor  using  calipers;  see  Figure  13).  This  finding 
demonstrates  selective  probe  bio-distribution  at  early  stages  of  disease  progression  and 
represents  a  significant  achievement  toward  the  goal  of  early  detection. 


Figure  13:  Tomographic  slice  (z  =  1.28  cm)  revealing  a 
1  mm  tumor  in  the  right  mammary  fat  pad  of  mouse. 


A2.5:  Evaluation  of  status  of  achieving  vear-2  goals 


As  previously  stated  we  have  already  made  significant  progress  toward  achieving  the  year  2 
goals.  We  are  confident  in  our  optimization  of  imaging  parameters  as  well  as  the  dose 
administration  of  the  NIR,  cathepsin-activatable  probe.  What  remains  ahead  for  year  2  is 
completion  of  the  screening  study.  To  continue  toward  completion  of  this  goal,  we  will  continue 
to  breed  mice  to  achieve  the  necessary  number  of  subjects.  As  stated  in  the  grant  proposal,  we 
will  screen  mice  as  they  develop  spontaneous  tumors  of  the  mammary  fat  pads.  We  fully 
anticipate  achieving  our  target  of  40  experimental  subjects  and  10  controls  ahead  of  schedule. 
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Additionally,  we  have  already  begun  correlative  work,  such  as  histological  examination  by 
hematoxylin  and  eosin  staining,  immunohistochemical  staining  for  cathepsins,  and  Western 
blotting  for  cathepsins,  to  confirm  the  presence  or  absence  of  disease  as  well  as  to  define  the 
sensitivity  thresholds  for  early  detection  by  FMT  and  reflectance  imaging  in-vivo.  Therefore,  we 
foresee  successful  completion  of  the  year-2  goals  ahead  of  schedule,  which  will  enable  us  to 
more  quickly  proceed  toward  accomplishing  the  goals  for  year  3. 

Aim  3  (Year  3):  To  translate  animal  study  findings  into  clinically -relev ant  detection  schemes. 

Consistent  with  the  time  plan  of  the  proposal  in  the  second  year  we  will  be  constructing 
phantoms  that  simulate  the  human  breast,  both  geometrically  and  in  terms  of  optical  properties 
and  we  will  performing  feasibility  studies  to  identify  clinical  utility. 

Key  Research  Accomplishments 


•  Implementation  of  a  system  that  combines  4  modes  of  illumination/detection  namely 
reflectance,  transillumination  tomosynthesis  and  tomography  for  use  in  years  2  &  3 

•  Construction  of  appropriate  phantoms  for  testing  imaging  performance. 

•  Establishment  of  appropriate  animal  breast  cancer  models. 

•  Establishment  of  sensitivity  limits  for  the  detection  of  fluorochromes  of  better  than 
lOnM  in  200pL  volumes  in  the  center  of  phantoms  with  breast  optical  properties. 

•  Establishment  that  tomographic  methods  are  superior  to  planar  imaging  methods  as  a 
function  of  depth,  optical  properties  and  heterogeneity  and  is  the  method  of  choice  for 
subsequent  studies  in  years  2  &  3.  Planar  imaging  methods  will  be  important  for 
scout  studies,  where  quantification  will  be  necessarily  obtained  from  tomography  or 
tomosynthesis. 

•  Development  of  algorithms  that  are  insensitive  to  background  fluorescence  and 
heterogeneity. 

•  In-vivo  imaging  of  the  first  animals  with  spontaneous  breast  cancer  using  protease 
activatable  probes  to  detect  cancer 

Reportable  Outcomes 

Journal  Publications  &  Proceedings: 

Ntziachristos  V,  Turner  G,  Dunham  J,  Windsor  S,  Soubret  A,  Ripoll  J  Shih  HA 

“Planar  fluorescence  imaging  using  normalized  data,”  Journal  of  Bio-optics,  in  press,  2005. 

Soubret  A,  Ripoll  J,  Ntziachristos  V,  “Accuracy  of  fluorescence  tomography  in  the  presence  of 

heterogeneities:  Study  of  the  normalized  Born  ratio,”  IEEE  Trans  Med  Imag,  in  press,  2005. 
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Soubret  A,  Windsor  SD,  Turner  GM,  Ripoll  J,  Ntziachiistos  V,  “Accuracy  of  fluorescence 
tomography  in  the  presence  of  background  fluorescence,”  in  preparation,  2005. 

Windsor  SD,  Shih  HA,  and  Ntziachiistos,  V,  “In-vivo  fluorescence  molecular  tomography  of 
mammaiy  adenocarcinomas  in  transgenic  mice  bearing  an  activated  c-neu  oncogene,”  U.S.  Army 
Era  of  Hope  Conference,  Philadelphia,  PA,  June  2005. 

Windsor  SD,  Shih  HA,  Weissleder  R,  Ntziachiistos  V,  “Simultaneous  imaging  of  protease 
expression  and  biodistribution  in  Her-2/neu  mice  using  dual -wavelength  fluorescence  molecular 
tomography,”  Society  of  Molecular  Imaging  Conference,  Cologne,  Germany,  Sept.  2005. 

Presentations: 

Windsor  SD,  “In-vivo  fluorescence  molecular  tomography  of  mammary  adenocarcinomas  in 
transgenic  mice  bearing  an  activated  c-neu  oncogene,”  U.S.  Army  Era  of  Hope  Conference, 
Philadelphia,  PA,  June  09,  2005. 

Conclusions 

In  this  year"  1  we  achieved  all  proposed  goals  and  met  all  milestones,  hi  addition  we  proceeded  in 
in-vivo  studies,  originally  planned  for  year  2  with  significant  success.  We  were  able  to  identify 
tumors  as  early  as  2-3mm  and  we  were  able  to  follow  their  growth  as  a  function  of  time.  In  years 
2  &  3  we  will  continue  studies  towards  understanding  the  targeting  ability  and  early  detection 
based  on  protease  probes  and  examine  the  potential  for  clinical  translation. 
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Accuracy  of  Fluorescent  Tomography 
in  thePresence  of  Heterogeneities: 
Study  of  theNormalized  Born  Ratio 

Antoine  Soubret*,  Jorge  Ripoll,  and  Vasilis  Ntziachristos 


Abstract — We  studied  the  performance  of  three-dimensional 
fluorescence  tomography  of  diffuse  media  in  the  presence  of 
heterogeneities.  Experimental  measurements  were  acquired  using 
an  imaging  system  consisting  of  a  parallel  plate-imaging  chamber 
and  a  lens  coupled  charge  coupled  device  camera,  which  enables 
conventional  planar  imaging  as  well  as  fluorescence  tomography. 
To  simulate  increasing  levels  of  background  heterogeneity,  we 
employed  phantoms  made  of  a  fluorescent  tube  surrounded 
by  several  absorbers  in  different  combinations  of  absorption 
distribution.  We  also  investigated  the  effect  of  low  absorbing 
thin  layers  (such  as  membranes).  We  show  that  the  normalized 
Born  approach  accurately  retrieves  the  position  and  shape  of  the 
fluorochrome  even  at  high  background  heterogeneity.  We  also 
demonstrate  that  the  quantification  is  relatively  insensitive  to  a 
varying  degree  of  heterogeneity  and  background  optical  proper¬ 
ties.  Findings  are  further  contrasted  to  images  obtained  with  the 
standard  Born  expansion  and  with  a  normalized  approach  that 
divides  the  fluorescent  field  with  excitation  measurements  through 
a  homogeneous  medium. 

Index  Terms — Absorption,  diffusion  theory,  fluorescence  tomog¬ 
raphy,  molecular  imaging. 


I.  Introduction 

FLUORESCENCE  tomography  of  diffuse  media  has  re¬ 
cently  materialized  as  a  quantitative  method  for  imaging 
novel  fluorescent  agents  with  functional  or  molecular  specificity 
[1],  Tomographic  approaches  have  been  developed  to  overcome 
significant  limitations  present  in  planar  optical  imaging  sys¬ 
tems,  i.e.,  systems  that  front  illuminate  the  animal  of  interest 
at  an  excitation  wavelength  with  a  broad  photon  beam  and 
record  fluorescent  images  of  the  animals  using  an  appropriate 
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high-pass  or  bandpass  filter.  Planar  systems  yield  images  that 
are  a  superposition  of  photons  emitted  from  different  depths 
and  scattered  by  the  different  layers  of  tissues.  In  this  config¬ 
uration,  depth  resolution  is  not  possible  and  quantification  of 
the  fluorescent  agent  is  further  complicated  by  the  nonlinear 
dependence  of  fluorescence  strength  on  depth  and  tissue  optical 
properties.  Conversely  tomography  utilizes  appropriate  photon 
propagation  models  that  account  for  nonlinear  effects  and 
multiple  photon  projections  and  can  independently  resolve  the 
three-dimensional  (3-D)  position  of  fluorochromes  and  their 
strength  and  concentration. 

Several  groups  have  offered  theoretical  and  experimental 
solutions  to  the  tomographic  problem  [2]— [18].  Inversion  is 
often  based  on  solutions  of  the  diffusion  equation  obtained 
either  analytically  or  numerically.  Tomographic  performance 
has  been  evaluated  with  simulated  data  [4],  [7],  [11],  [12],  and 
[19]— [23]  and  with  experimental  measurements  from  phantoms 
[3],  [5],  [17],  [24].  Most  algorithmic  evaluations,  however,  have 
been  performed  with  relatively  homogeneous  media;  typically 
a  small  number  of  fluorochromes  is  immersed  in  an  otherwise 
homogeneous  fluid  with  average  optical  properties  similar  to 
those  of  a  target  tissue  (for  example  the  human  breast  or  a  small 
animal).  Therefore,  a  common  criticism  of  the  technique  relates 
to  its  accuracy  and  performance  in  media  with  high  background 
optical  heterogeneity  such  as  tissues. 

Developing  robust  methods  that  work  well  at  highly  hetero¬ 
geneous  diffuse  media  is  important  for  in  vivo  imaging.  We  have 
previously  reported  on  a  particular  methodology,  termed  the 
normalized  Bom  ratio  that  demonstrated  good  imaging  perfor¬ 
mance  in  vivo.  The  corresponding  imaging  technology,  termed 
fluorescence  molecular  tomography  (FMT)  [25],  is  based  on 
using  self  calibrated  data  that  divide  the  fluorescence  measure¬ 
ments  with  corresponding  measurements  at  the  excitation  wave¬ 
length  [17].  This  normalization  scheme  can  work  with  analytical 
or  numerical  forward  solvers  and  offers  significant  experimental 
advantages  since  it  is  independent  of  source  strengths,  detector 
gains  and  coupling  efficiency  to  tissue  [17].  Furthermore,  it  min¬ 
imizes  sensitivity  to  theoretical  errors  such  as  modeling  of  the 
boundary  conditions  or  knowledge  of  the  exact  tissue  optical 
properties.  This  normalization  scheme  was  pivotal  to  the  devel¬ 
opment  of  the  first  in  vivo  fluorescence  observations  of  protease 
up-regulation  deep  in  mouse  brains  [26]  and  of  tumor  treatment 
in  intact  animals  [27],  using  appropriate  fluorescent  probes. 

The  normalized  Bom  ratio  has  been  shown  insensitive  to  a 
small  number  of  absorbers  in  the  past  [17]  using  experimental 
measurements.  Furthermore,  Roy  et  al.  [19]  demonstrated  that 
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Fig.  1.  Experimental  setup. 


high  levels  of  noise  could  be  added  to  simulated  data  without 
drastically  changing  the  reconstruction  performance  when  using 
this  normalization  method.  Similarly,  Eppstein  et  al.  [23]  have 
shown  relative  insensitivity  of  fluorescence  reconstructions  to 
the  spatial  variation  of  optical  properties  using  simulated  data. 

The  work  herein  focused  on  experimentally  exploring  the  per¬ 
formance  of  the  normalized  Born  ratio  with  increasing  back¬ 
ground  absorption  heterogeneity.  We  employed  phantoms  con¬ 
taining  a  varying  number  of  diffuse  absorbing  tubes  with  dif¬ 
ferent  absorption  coefficients  surrounding  a  fluorescent  diffuse 
tube.  While  it  was  hypothesized  that  the  normalized  Born  ratio 
could  perform  well  at  increasing  background  absorption,  our  re¬ 
sults  demonstrated  an  unexpected  robustness  even  at  very  high 
degree  of  heterogeneity,  despite  the  linear  inversion  approach 
used.  To  contrast  these  findings,  we  also  examined  the  perfor¬ 
mance  of  1)  the  standard  Bom  reconstruction  scheme  and  2)  a 
normalization  scheme  that  uses  calibrated  data  based  on  mea¬ 
surements  through  a  homogeneous  medium  in  the  absence  of 
the  absorbing  tubes.  Section  II,  describes  the  experimental  setup 
and  the  different  heterogeneities  considered.  Section  III  presents 
the  forward  model  and  the  inversion  algorithm  employed.  Sec¬ 
tion  IV,  describes  our  findings  and  finally  Section  V  discusses 
the  results  and  their  implications  to  in  vivo  imaging. 

II.  Experimental  Setup  and  Phantoms 
A.  Experimental  Setup 

A  schematic  of  the  FMT  acquisition  system  is  shown  in  Fig.  1 
and  a  detailed  description  can  be  found  in  [16].  A  672  nm  wave¬ 
length  constant  wave  (CW)  laser  is  routed  into  a  two  channel 
optical  switch.  This  switch  allows  the  system  to  be  used  either 
as  a  reflectance  or  tomographic  imager.  In  reflectance  mode,  the 
laser  beam  was  expanded  with  a  system  of  lenses  and  front-il- 
luminates  the  object  placed  in  the  imaging  chamber  though  a 
window  with  anti-reflection  coating  for  the  near  infrared  re¬ 
sponse,  In  tomographic  mode,  a  programmable  optical  switch 


sequentially  directed  light  in  each  of  the  46  fibers  arranged  in 
a  space-filling  grid  pattern  used  for  tomography  purposes.  The 
fibers  were  arranged  on  a  moving  plate  that  can  establish  a  slab- 
geometry  of  different  thicknesses;  here  in  1.7  cm.  The  chamber 
can  be  filled  with  an  optical  matching  fluid  with  absorbing  and 
scattering  properties  similar  to  that  of  tissues.  The  light  trans¬ 
mitted  through  the  diffuse  medium  is  then  recorded  with  the 
charge  coupled  device  (CCD)  camera  at  both  the  excitation  and 
fluorescent  wavelength  using  appropriate  filters  (for  details  see 
[16]). 

B.  Phantom  Geometries 

Different  configurations  of  background  heterogeneity  were 
constructed  by  surrounding  a  single  fluorochrome  with  multiple 
absorbers  of  increasing  density  and  absorption.  The  fluorescent 
tube  was  made  of  a  3-mm-diameter  translucent  plastic  tube  con¬ 
taining  1  pM  of  Cy  5.5  in  a  solution  of  1%  intralipid  and  50  ppm 
of  India  ink.  Absorbers  were  made  from  the  same  tube  filled 
with  1%  intralipid  and  different  amounts  of  India  ink.  All  tubes 
were  placed  in  position  using  a  solid  frame  that  does  not  inter¬ 
fere  with  the  field  of  view  and  immersed  in  a  solution  of  1  % 
intralipid  and  50  ppm  of  India  ink  yielding  optical  properties 
pPa  —  0.4  cm-1  and  p's  =  8  cm-1  which  are  characteristic  of 
small  animals.  Therefore,  the  fluorescent  tube  and  background 
medium  had  the  same  optical  properties. 

Three  distinct  sets  of  experiments  were  performed.  The  first 
set  a)  consisted  of  12  configurations  al-al2  as  shown  in  Fig.  2. 
The  number  of  plastic  tubes  in  the  chamber  is  kept  constant  to 
eight  except  for  the  configuration  al2  and  the  attenuation  of  the 
measured  signal  from  one  configuration  to  another  is  only  due 
to  the  different  concentration  of  ink  in  the  tubes.  The  second 
set  of  experiments  b)  consisted  of  seven  configurations  and  it 
is  similar  to  the  first  set  of  experiments  al-a7  with  the  excep¬ 
tion  that  tubes  with  the  same  optical  properties  with  the  back¬ 
ground  medium  are  removed  from  the  medium.  Therefore,  the 
effect  of  the  plastic  tube  heterogeneity  can  be  assessed.  Those 
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Fig.  2.  Horizontal  slices  (xz  plane)  of  the  different  phantoms  made  of  a  fluorochrome  tube  surrounded  by  absorbing  tubes.  Three  sets  of  experiments  a),  b),  and 
c)  have  been  used  to  perform  the  experiments.  The  3-D  view  of  the  phantoms  al-al  1  and  b7  is  shown  on  the  view  LI).  The  legend  indicating  the  concentration  of 
ink  in  the  different  tubes  is  depicted  on  the  view  L2). 


seven  configurations  are  numbered  bl-b7  in  Fig.  2.  The  third 
set  of  experiments  c)  imaged  a  fluorescent  tube  placed  super¬ 
ficially  in  the  diffuse  fluid,  virtually  in  contact  with  the  inner 
side  of  the  imaging  window.  Two  configurations  were  examined, 
one  for  imaging  the  tube  placed  in  a  homogenous  background 
medium  and  the  other  with  the  tube  surrounded  by  seven  ab¬ 
sorbers.  These  two  configurations  are  named  cl  and  c7  in  Fig.  2. 

Fig.  3  shows  an  example  of  the  phantoms  constructed;  in 
particular  it  shows  reflectance  images  (photographs)  of  the 
phantom  set  bl  and  set  b7  of  Fig.  2  when  there  is  no  back¬ 
ground  matching  fluid  (ink-intralipid  solution)  in  the  chamber. 
Photographs  at  the  excitation  and  the  fluorescent  channels  are 
shown  respectively.  As  seen,  almost  no  fluorescent  signal  is 
detectable  in  the  configuration  b7  as  compared  to  the  config¬ 
uration  bl  where  no  absorbers  are  present.  When  filling  the 
chamber  with  the  diffuse  medium  the  fluorescent  tubes  cannot 
be  directly  observed  with  reflectance  imaging  due  to  high 
diffusion  and  can  be  revealed  only  tomographically. 

III.  Data  Normalization  and  Reconstruction 
Algorithm 

A.  Forward  Model 

The  propagation  of  near-infrared  light  in  tissue  is  modeled 
by  the  diffusion  equation  which  has  been  shown  to  work  well 


when  sources  and  detectors  are  separated  by  more  than  a  few 
mean  free  path  lengths  [28].  In  our  experiment,  more  than  13 
mean  free  path  separate  the  sources  and  the  window  imaged 
by  the  CCD  camera.  When  the  fluorescent  tube  is  inserted  in 
the  chamber,  the  fluorescent  field  Uftao(rd,  rs)  detected  by 
the  CCD  camera  at  the  exit  plane  is  related  to  the  absorption 
of  the  fluorochrome  HaUO(r l)  by  0),  shown  at  the  bottom  of 
the  page,  where  n{r{)  =  t?/Ual'0(ri)>  V  the  quantum  yield  of 
the  fluorochrome,  rs,  i\i  are  the  location  of  the  source  and  the 
detector,  ©3rc(?Y)  is  the  intensity  of  the  source  at  the  fiber  exit, 
©juo  the  attenuation  of  the  fluorescent  filter,  ©det(rd)  accounts 
for  the  detectors  gains,  Qglut  is  the  detector  quantum  efficiency 
at  the  fluorescence  wavelength,  v  is  the  speed  of  light  in  the 
medium  and  DAflu°  =  v/3  p/s  the  diffusion  coefficient  at  the 
fluorescence  wavelength.  The  Green’s  function  Gqc  (?’i ,  rs) 
describes  the  propagation  of  the  light  from  the  source  to  the 
fluorochrome  and  GQfluo(rt;,  ri)  from  the  fluorochrome  to  the 
detector.  These  functions  can  be  calculated  assuming  a  homo¬ 
geneous  background  medium  or  iteratively  using  nonlinear 
inverse  algorithms  based  on  finite-differences  or  finite-element 
methods  [7],  [8],  [23],  [29].  In  this  paper,  we  demonstrated 
imaging  performance  using  a  first-order  perturbative  solution 
of  the  diffusion  equation  assuming  a  homogeneous  background 
[1 1],  Boundary  conditions  were  implemented  using  the  extrap¬ 
olated  boundary  conditions  [30]— [32]  and  the  method  of  image 
sources  [33]  with  ten  source  pairs. 


rf3ri0det(rd)G^(rrf!n)- 


•  n(ri  )Gq' axc  (r  i ,  rs )  ©src  (rs ) 
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Fig.  3.  Reflectance  images  when  only  the  fluorochrome  tube  is  present  in  the  chamber  (configuration  bl)  and  when  the  fluorochrome  tube  is  surrounded  by  7 
absorbers  tubes(configuration  b7).  Images  are  for  the  excitation  channel  and  the  fluorescent  channel.  The  units  on  the  axis  x  and  y  are  in  centimeters. 


B.  Data  Normalization 


where 


The  use  of  (1)  is  experimentally  complicated  since  the  pa¬ 
rameters  Qg1"0,  0det (rd),  ©si -c{rs)  have  to  be  explicitly  deter¬ 
mined.  In  addition,  the  measurements  can  be  complicated  by 
unknown  coupling  factors  of  the  optical  system  onto  tissue.  We 
have  previously  suggested  the  normalized  Bom  ratio,  which  di¬ 
vides  fluorescence  measurements  with  corresponding  measure¬ 
ments  at  the  excitation  wavelength  [17],  in  order  to  simplify 
experimental  issues.  This  is  because  the  normalized  Bom  ratio 
is  effectively  independent  of  the  individual  gains  of  source  de¬ 
tector  pairs  and  coupling  coefficients.  Herein,  we  examine  its 
performance  with  diffusive  media  of  heterogeneous  absorption 
distribution. 

To  write  the  normalized  Born  ratio,  we  first  decompose  the 
excitation  field  using  Green  functions  as 


UeKC(rd,rs)  =  G'oe"c(rd,rs)  0src (rs) 

(2) 

where  Q£ *'  is  the  detector  quantum  efficiency  at  the  excitation 
wavelength  and  0^c  is  the  attenuation  of  the  filter  at  the  exci¬ 
tation  wavelength.  Then,  the  normalized  Bom  ratio  U^(rd,  rs ) 
is  obtained  by  dividing  (1)  by  (2),  i.e.. 


U<l\rd,r.)=  VaUrj^s) 

Uexcy  dt^s) 


/d3riGofIuo(r</,ri)n(ri)Go“c(ri,r1,) 

Go**c{rdirs) 


Qgf,uo  0f 0  « 

gAexc  Qexc  jjX £i, 


The  single  calibration  parameter  ao  depends  on  generic  system 
characteristics  (i.e.,  filters’  and  fibers’  attenuation,  laser  power, 
etc.)  and  can  be  simply  determined  experimentally  by  measure¬ 
ment  of  a  fluorochrome  with  known  concentration. 

In  some  cases  imperfect  isolation  of  the  fluorescent  filter  may 
result  in  small  remnant  signals  ©bleed  of  the  excitation  signal 
propagating  in  the  fluorescent  channel.  In  such  occurrence,  this 
“bleed-through”  signal  is  subtracted  from  the  fluorescence  field, 
i.e., 


[' Um{rd,ra ) 


% 


measured 

)  ©bleed  U excis'd :  Fs ) 


Gexc  (Vd ;  ) 


(4) 


J  measured 


The  coefficient  ©bleed  can  be  experimentally  determined  by  il¬ 
luminating  the  chamber  filled  with  any  nonfluorescent  diffusing 
medium  and  calculating  the  ratio  of  counts  measured  at  the  fluo¬ 
rescent  channel  over  the  counts  measured  at  the  intrinsic  channel 
[16]. 

As  it  will  appear  in  the  following,  this  Born  ratio  does  not 
only  offer  convenient  calibration  but  also  significantly  mini¬ 
mizes  the  sensitivity  to  the  presence  of  background  optical  het¬ 
erogeneities.  The  following  three  different  data  schemes  were 
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examined  here  in  relation  to  their  performance  with  increased 
background  heterogeneity. 

•  Normalization  I,  which  employed  the  normalized  Born 
ratio  as  shown  in  (4).  The  excitation  light  field  used  was 
obtained  from  the  heterogeneous  medium  containing  the 
fluorochrome  and  the  different  combinations  of  absorbers. 

•  Normalization  II,  which  employed  a  normalized  Born 
field  U^ll\rd,  rs),  whereas  the  excitation  light  field 
U£r°s\rd,  rs )  is  measured  from  a  homogeneous  back¬ 
ground  medium,  i.e.,  from  the  chamber  filled  only  with  a 
homogeneous  solution  of  intralipid  and  ink  in  the  absence 
of  absorbers  or  fluorochromes.  This  fluid  was  identical  to 
the  fluid  used  as  a  background  medium  in  the  heteroge¬ 
neous  case.  The  field  employed  in  this  normalization  can 
be  written  similarly  to  (4)  as 


'uM{rd,r.j 


measured 


EWrrf,r,)  -  &^utmos)(rd,rs 

Utm°s\rd,rs) 


(5) 


J  measured 


•  Normalization  III,  which  employed  a  standard  Bom  field 
[f(niHrd,rs)  (see  [11]).  This  field  is  normalized  only 
with  source  strength  <?o(rs)  at  each  position  rs  since  the 
use  of  individual  fibers  in  the  system  described  above 
often  leads  to  variability  in  optical  transmittance  from 
fiber  to  fiber  due  to  differences  in  coupling  losses  in  the 
optical  switch  and  elsewhere.  The  field  employed  in  this 
normalization  can  be  written  as 


”  J  measured 

tfruofa, rs)  -  eUeedUtm°s)(rd,rs 
9o(rs) 


(6) 


measured 


The  factor  go  ( ra )  is  calculated  as  the  maximum  average 
intensity  of  photon  fields  propagating  though  the  homo¬ 
geneous  medium  for  the  different  sources  employed  so 
that  go(rs)  =  ma,xrd(Ui^°m°s) (rd, rs)).  This  normal¬ 
ization  works  on  the  basis  that  go{i's)  is  proportional  to 
the  intensity  of  the  source  @SVc(rs)  located  at  rs.  This 
method  normalizes,  therefore,  the  measurements  to  the 
source  strength  but  it  does  not  take  into  account  the  spatial 
dependence  of  the  measured  excitation  field  since 
only  depends  on  the  source  position  rs  and  not  on  the  lo¬ 
cation  rd  of  the  detector. 


Then  the  discretized  linear  forward  problem  can  be  written  as 
[^measured  =  «0  ]T  Wi  fni ■  (8) 

3 

As  we  do  not  take  into  account  the  heterogeneities  in  the  eval¬ 
uation  of  the  Green  functions  the  same  weight  matrix  is  used 
when  the  Normalization  II  is  considered 

[tf  ^measured  =  «0  ^  Wifni  ■  (9) 

3 

For  implementing  the  normalization  III,  fluorescent  data  are 
normalized  only  with  the  corresponding  source  strengths.  Con¬ 
sequently  the  weight  function  is  not  normalized  in  this  case  by 
the  source-detector  Greens  function  and  the  forward  problem  is 
written  as 


W |  =  6VGx0^(r'dy)Gt'(r3,rl) 


and  the  corresponding  forward  problem  is  written  as 

[^III)]meaSUred  =  ao^^nj 

j 


where  «o  is  the  corresponding  normalization  factor  for  the  Bom 
normalization 


d0  =  Qe  © 


•Wo  C)fluO  V 


E  '■'f  £>W o  ' 


The  3-D  mesh  used  for  the  reconstructions  was  Nx  =  29, 
Ny  —  27,  Nz  =  18  for  a  field  of  view  of  2.8  cm  X  2.6  cm 
with  a  depth  equal  to  the  thickness  of  the  slab  (i.e.,  1.7  cm). 
The  voxel  size  was  0.1  cm  x  0.1  cm  x  0.1  cm.  A  subsampled 
set  of  measurements  was  extracted  from  the  512  x  512  pixels 
image  acquired  by  the  CCD  camera  by  averaging  square  re¬ 
gions.  This  process  results  in  a  total  of  30  x  30  detectors  avail¬ 
able  for  each  measurement.  For  each  source  detector  pair  (which 
consists  of  30  x  30  X  46  =  41  400  values),  a  threshold  was  ap¬ 
plied  to  select  only  the  measurements  with  a  high  signal  to  noise 
ratio  (i.e.,  around  10000-15  000  values  are  effectively  used). 
The  fluorescence  rij  is  obtained  in  calculating  the  Moore-Pen- 
rose  pseudoinverse  with  an  algebraic  reconstruction  technique 
(ART)  with  a  random  projection  access  order  [34],  The  number 
of  ART  iterations  was  fixed  to  20.  Several  relaxation  parame¬ 
ters  have  been  tested,  and  as  noticed  by  Herman  [35]  a  strong 
under-relaxation  scheme  gives  the  best  result  in  terms  of  solu¬ 
tion  stability.  In  the  following  reconstructions,  we  have  used  the 
value  0.03  for  the  relaxation  parameter. 


C.  Reconstruction 


To  solve  the  inverse  problem,  we  discretized  the  integral 
equation  (3)  and  evaluated  the  fluorescent  term  n(ri)  on  a  3-D 
grid  of  size  Nx,  Ny,  Nz.  For  each  source-detector  pair  (rd,  rls) 
and  each  point  (r,)j=i...jvI  jvyAt.  on  the  grid,  the  weight  func¬ 
tion  W-jB  connecting  the  measurement  U-13  =  UnB(rd,r\) 
to  the  unknowns  nj  =  n{r^)  is 


W?3  =  8V 


(7) 


where  S V  is  the  volume  of  each  voxel  of  the  grid. 


D.  Image  Analysis 

For  quantification  purposes  the  reconstructed  values  reported 
for  the  fluorescent  tube  were  calculated  as  a  summation  of  the 
fluorescence  intensities  reconstructed,  above  a  threshold,  for  all 
voxels  contained  within  a  selected  volume  of  interest  (VOI).  The 
threshold  was  empirically  set  to  5%  of  the  maximum  value  re¬ 
constructed  for  the  configuration  al,  for  each  of  the  three  nor¬ 
malizations,  and  was  kept  constant  for  subsequent  reconstruc¬ 
tions.  The  VOI  selected  contained  all  voxels  reconstructed,  ex¬ 
pect  for  the  ones  belonging  to  the  first  two  and  the  last  two  ver¬ 
tical  (coronal)  slices,  i.e  the  slices  that  are  parallel  to  the  x-y 
plane  and  are  respectively  close  to  the  source  plane  and  the 
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Fig.  4.  Vertical  (coronal)  slices  along  the  x-y  plane  at  the  depth  z  =  0.65  cm  showing  the  reconstructed  fluorochrome  tube  for  the  configurations  set  al,  a2,  a4,  a7, 
a8,  all,  al2.  The  first  row  I)  has  been  obtained  with  the  Normalization  I,  the  second  row  II)  with  the  Normalization  II  and  the  third  row  III)  with  the  Normalization 
III  as  described  in  Section  III-A.  For  each  row,  the  fluorescence  reconstructed  values  for  the  configuration  al-al2  have  been  divided  by  the  maximum  value  for 
the  configuration  al. 


imaging  window.  This  is  because  some  artifacts  tend  to  appear 
close  to  the  sources  and  detectors.  These  artifacts  were  thus  ex¬ 
cluded  from  the  quantification  calculations. 

To  evaluate  quantification  errors  as  a  function  of  normaliza¬ 
tion  method  and  configuration,  one  of  the  configurations  exam¬ 
ined  was  selected  as  a  reference  and  the  corresponding  summed 
fluorescence  intensity  reconstructed  was  fixed  to  1  (in  arbitrary 
units).  All  other  reconstructed  summed  intensities  were  divided 
with  this  reference  value,  to  record  the  relative  change  in  quan¬ 
tification,  as  result  of  the  normalization  method  selected  and 
the  absorber  configuration  examined.  Results  are  plotted  by  se¬ 
lecting  two  different  configurations  as  references,  i.e.,  the  con¬ 
figuration  al,  a6  for  the  set  of  experiment  a)  and  bl,  b4  for  the 
set  b). 

To  assess  overall  image  quality,  each  reconstructed  image 
was  attributed  with  a  score  (from  1  to  10)  based  on  localiza¬ 
tion  and  shape  accuracy.  These  parameters  were  visually  eval¬ 
uated  from  two  imaging  scientists  at  the  post-graduate  level 
trained  for  more  than  five  years  to  assess  imaging  quality  of 
diffusive  images  by  superimposing  reconstructed  images  on  the 
tube  photographs  obtained  under  identical  positional  parameters 
and  with  the  same  optical  system  (prior  to  adding  surrounding 
absorbers  and  the  diffuse  fluid).  This  scoring  method  was  used 
for  evaluating  the  relative  qualitative  imaging  performance  of 
the  different  normalization  methods  employed. 

IV.  Experimental  Results  and  Reconstruction 
A.  Imaging  at  the  Presence  of  Absorbing  Heterogeneity 

Fig.  4  shows  vertical  reconstructed  slices  (coronal  slices)  at 
the  depth  z  =  0.65  cm  of  the  3-D  reconstructed  tube  for  the 
configurations  al,  a2,  a4,  a7,  a8,  all,  and  al2.  Each  row  in  the 
image  corresponds  to  reconstructed  data  obtained  with  the  Nor¬ 
malizations  I,  II,  and  III,  respectively.  The  results  shown  for  dif¬ 
ferent  configurations  are  scaled  using  the  configuration  al  as 


reference  for  each  normalization  scheme  I,  II,  and  III.  Calibra¬ 
tion  of  the  constants  Qq  and  &o  can  be  used  to  convert  the  data  to 
concentration  units;  however,  this  was  not  practiced  herein  since 
it  does  not  bear  relevance  to  the  particular  goals  of  the  study. 

Evident  on  the  images  is  that  the  fluorescent  tube  is  well  re¬ 
constructed  at  all  levels  of  absorption  variation  attempted  when 
employing  the  normalization  I.  The  tube  position  and  size  are 
well  reconstructed  as  background  absorption  heterogeneity  is 
increased.  The  reconstructions  produced  using  the  normaliza¬ 
tion  II  are  also  fairly  accurate  but  the  shape  of  the  tube  is  dis¬ 
torted  and  right  shifted  for  high  absorption  as  clearly  visible  on 
the  set  al2.  Finally,  imaging  with  the  normalization  III  attains 
the  largest  distortion  since  the  reconstructed  images  are  strongly 
affected  in  shape,  quantification  and  position.  In  particular,  the 
fluorescence  intensity  reconstructed  is  stronger  at  the  center  of 
the  reconstructed  field.  Furthermore,  the  asymmetry  of  the  con¬ 
figuration  a  12,  where  the  absorbers  are  located  only  on  one  side 
of  the  fluorochrome,  is  evident  on  the  corresponding  reconstruc¬ 
tion. 

All  schemes  examined  yielded  variation  on  the  fluorescence 
intensity  reconstructed.  Fig.  5  summarizes  the  quantification 
variation  for  the  different  configurations  examined  and  the 
data  normalization  employed.  When  configuration  al  was  used 
as  reference,  the  subsequent  reconstructed  data  demonstrated 
a  gradual  overestimation  of  the  fluorescence  intensity  recon¬ 
structed  for  normalization  I.  This  overestimation  is  of  the  order 
of  10%  for  the  first  seven  cases  (set  al  to  a7).  As  the  absorption 
further  increases,  the  reconstructed  value  is  overestimated  by 
up  to  15%.  In  contrast,  an  underestimation  is  observed  for 
normalizations  II  and  III.  Here,  the  underestimation  exceeds 
60%  for  the  most  absorbing  cases. 

Fig.  5(b)  re-plots  this  data,  assuming  configuration  a6  as 
a  reference.  This  plot  represents  better  the  in  vivo  situation 
where  some  degree  of  heterogeneity  is  always  present.  Based 
on  these  plots  it  is  noted  that  Normalization  I  generally  leads 
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Fig.  5.  Relative  change  in  fluorescence  strength  reconstructed  for  the  different  data  normalizations  as  a  function  of  configurations  al-al2.  (a)  Relative 
quantification  change  as  a  function  of  background  heterogeneity  assuming  configuration  al  as  reference,  (b)  Re-plot  of  the  data  in  (a)  assuming  the  configuration 
a6  as  reference. 
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Fig.  6.  Vertical  (coronal)  slices  along  the  x-y  plane  at  the  depth  z  =  0.65  cm  and  horizontal  (axial)  slices  along  the  x-z  at  the  position  y  =  0  cm  showing  the 
reconstructed  fluorescent  tube  for  the  configurations  bl-b3.  The  first  row  I)  has  been  obtained  with  the  normalization  I,  the  second  row  II)  with  the  normalization 
II  and  the  row  III)  with  the  normalization  III.  For  each  row,  the  images  have  been  divided  by  the  maximum  value  for  the  configuration  bl. 


to  a  fluorochrome  strength  over-estimation  for  increasing  back¬ 
ground  absorption  heterogeneity  and  to  an  underestimation  for 
decreasing  background  heterogeneity.  Correspondingly,  Nor¬ 
malizations  II  and  III  result  in  overestimation  for  decreasing 
background  absorption  heterogeneity  and  to  an  underestimation 
for  increasing  background  heterogeneity.  This  result  points  to 
the  practical  and  intuitive  conclusion  that  algorithm  calibra¬ 
tions  should  be  performed  in  heterogeneous  media  to  reduce 
quantification  errors  for  background  absorption  variations.  We 
have  also  performed  reconstructions  using  a  forward  model  that 
assumed  different  absorption  coefficients  that  match  the  av¬ 
erage  optical  properties  for  each  configuration  examined.  This 
yielded  similar  observations,  especially  for  the  normalization 
methods  I  and  II  due  to  the  relative  insensitivity  of  the  ratio 
methods  to  the  exact  knowledge  of  the  optical  properties. 

B.  Effect  of  Semi-Translucent  Inteifaces 

Fig.  6  shows  the  vertical  reconstructed  slices  (coronal  slices) 
at  the  depth  z  =  0.65  cm  and  horizontal  slices  (axial  slices)  at 
the  position  y  =  0  cm  for  sets  bl  -  b3.  Similar  imaging  per¬ 
formance  was  obtained  for  the  other  configurations  as  well.  All 


images  are  referenced  to  the  configuration  bl,  i.e.,  each  of  the 
image  shown  was  divided  by  the  maximum  value  obtained  for 
the  configuration  bl.  Similarly  to  Fig.  4,  the  results  depict  that 
the  Normalization  I  reconstruction  accurately  reconstructs  the 
location  and  shape  of  the  tube  on  the  coronal  views.  Normaliza¬ 
tion  II  and  III  are  more  sensitive  to  the  background  optical  vari¬ 
ation  since  the  yield  less  accurate  shape  reconstructions  com¬ 
pared  to  the  ones  obtained  with  Normalization  I.  The  effects  of 
the  background  heterogeneity  seen  herein  are  more  significant 
than  the  ones  seen  in  Fig.  4  since  measurements  are  affected  by 
both  absorption  and  refraction  index  changes  from  one  config¬ 
uration  to  another  due  to  the  introduction  of  the  semi-translu- 
cent  plastic  tubes.  Interestingly,  the  reconstruction  shows  less 
distortion  for  higher  degree  of  heterogeneity  for  all  normaliza¬ 
tion  schemes  (results  not  shown)  possibly  because  of  the  relative 
symmetry  that  is  restored  in  the  medium  when  more  absorbing 
tubes  are  added. 

These  effects  are  also  evident  on  the  corresponding  axial  im¬ 
ages  (XZ  cross-sections)  shown  in  Fig.  6.  Since  the  scanner  op¬ 
erates  on  a  slab  geometry,  the  resolution  along  z  is  compromised 
in  relation  to  the  resolution  achieved  in  the  x-y  plane.  There- 
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Fig.  7.  Relative  change  in  fluorescence  strength  reconstructed  for  the  three  normalization  schemes  studied  as  a  function  of  configurations  bl-b7.  (a)  Results 
plotted  using  the  configuration  bl  as  reference,  (b)  Results  plotted  using  configuration  b4  as  the  reference. 


fore,  the  tube  is  reconstructed  elongated  using  Normalization  I 
but  its  x-y  location  is  well  resolved.  However,  Normalization  II 
and  HI  give  notable  distortion  in  the  tube  reconstructed. 

Fig.  7  depicts  the  quantification  trends  observed  for  the 
different  heterogeneity  configurations  examined.  Two  configu¬ 
rations  (bl  and  b4)  were  used  as  reference,  i.e.,  the  fluorescence 
strength  calculated  for  bl  or  b4  was  used  to  divide  all  other 
quantification  calculation  obtained.  Findings  indicate  a  similar 
performance  as  the  one  seen  in  Fig.  5.  The  normalization  I 
over-estimates  the  fluorescence  strength  for  increasing  back¬ 
ground  heterogeneity  and  underestimates  it  for  decreasing 
heterogeneity.  The  opposite  is  observed  for  Normalizations  II 
and  III.  However,  the  errors  in  quantification  are  stronger  than 
in  the  cases  examined  in  Fig.  5  due  to  the  added  effect  of  the 
semi-translucent  interfaces.  In  particular  when  the  reference 
configuration  bl  is  used,  we  observe  that  the  normalization 
I  yielded  overestimation  by  more  than  35%  when  more  than 
four  absorbers  were  present.  Correspondingly,  we  observed 
an  underestimation  of  up  to  65%  compared  to  the  reference 
reconstruction  when  using  normalization  II  and  III. 

C.  Effect  of  Depth 

To  examine  the  dependence  of  the  methods  to  depth,  we 
placed  the  tube  adjacent  to  the  front  imaging  window  and 
examined  the  reconstruction  performance  assuming  a  homo¬ 
geneous  background  and  a  highly  heterogeneous  background 
employing  seven  objects  as  shown  in  Fig.  2(c).  The  reconstruc¬ 
tions  using  the  normalization  I  are  shown  in  Fig.  8.  In  this  case, 
the  background  heterogeneity  bears  insignificant  effects  on  the 
image  quality  both  in  terms  of  shape  and  quantification. 

V.  Discussion  and  Conclusion 

In  this  paper,  we  examined  the  effect  of  background  optical 
heterogeneity  on  fluorescence  tomography  performance.  This 
work  is  strongly  linked  to  imaging  of  fluorochromes  in  tissues 


No  absorbers  7  absorbers 


x 


Fig.  8.  Fluorescence  reconstruction  using  the  normalization  I  when  the 
fiuorochrome  is  placed  superficially,  i.e.,  close  to  the  detector  plane  (see  Fig.  2). 
The  first  column  is  the  reconstruction  for  the  configuration  cl  without  any 
absorbers  and  the  second  column  for  the  configuration  c7  when  7  absorbers  are 
present.  The  vertical  (coronal)  slice  along  the  x-y  plane  is  at  depth  z  =  1 .6  cm. 
The  reconstructed  values  are  in  arbitrary  units. 

and  comes  to  address  concerns  associated  with  the  in  vivo  per¬ 
formance  of  fluorescence  tomography,  in  particular  when  using 
linear  solutions  derived  on  the  assumption  of  photon  propaga¬ 
tion  in  homogenous  diffuse  media. 

The  results  demonstrate  robust  performance  even  at  high 
background  optical  variations,  especially  when  employing 
the  Normalization  I.  This  robustness  was  demonstrated  using 
absorbing  variations  over  a  wide  variety  of  optical  properties 
and  spatial  distributions.  Good  imaging  performance  was  also 
demonstrated  for  the  case  where  the  effect  of  semi-translucent 
plastic  tubes  was  progressively  introduced  in  the  measurements. 
The  normalization  II  demonstrated  good  imaging  quality  only 
when  the  background  absorption  was  varied,  but  it  was  more 
sensitive  to  the  changes  introduced  by  the  translucent  tubes 
compared  to  the  performance  of  normalization  I.  Generally, 
these  normalized  methods  are  quite  insensitive  to  the  back¬ 
ground  optical  properties  used  in  the  reconstructions.  It  is, 
however,  possible  that  these  results  could  be  further  improved 
by  adapting  the  optical  properties  used  in  the  reconstruction  to 
the  exact  “average”  optical  properties  of  the  tissue  measured.  In 
addition,  the  results  obtained  using  the  normalization  II  could 
be  improved  if  a  separate  matching  fluid  was  employed,  with 
the  average  optical  properties  of  each  of  the  different  configu¬ 
rations  examined.  However,  this  practice  would  be  tedious  to 
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implement  in  in  vivo  measurements.  Finally,  the  normalization 
III  gave  overall  the  worse  imaging  performance  in  terms  of 
localization,  shape  accuracy,  and  reconstruction  owing  to  the 
spatial  variation  of  the  background  optical  properties  and  other 
experimental  and  theoretical  uncertainties. 

Our  observations  are  based  on  CW  measurements  due  to  their 
relevance  to  a  scanner  developed  for  small  animal  FMT  imaging 
[16].  Absorption  heterogeneities  cannot  be  distinguished  with 
high  contrast  from  scattering  heterogeneities  [36],  [37]  when 
using  CW  measurements  purely;  therefore,  the  study  of  absorp¬ 
tion  heterogeneities  in  this  paper  should  map  the  effects  of  con¬ 
trast  from  attenuation  in  general.  For  this  reason  instead  of  re¬ 
peating  measurements  with  scattering  diffuse  inhomogeneities 
we  opted  to  examine  the  effects  of  the  virtually  semi-diffusive 
plastic  surfaces  composing  the  tubes  employed  in  the  study.  This 
heterogeneity  can  be  seen  as  simulating  different  tissue  mem¬ 
branes,  layers  of  adipose  tissue  or  even  fluid  layers  as  common 
in  brain  measurements  due  to  the  cerebrospinal  fluid  [38],  [39] 
although  the  relative  geometry  of  these  semi-diffusive  surfaces 
will  also  influence  the  result.  It  was  experimentally  observed 
that  the  effects  of  the  tubes  yielded  significant  changes  in  the 
reconstructed  data,  and  these  effects  were  more  evident  in  the 
reconstructions  using  the  normalization  II  and  III.  However,  it 
is  expected  that  most  in  vivo  measurements  do  not  have  such 
high  density  of  semi-diffusive  regions,  especially  in  comparison 
to  the  highly  heterogeneous  configurations  (b5-b7)  examined 
herein 

The  normalization  I  demonstrated  better  quantification 
accuracy  than  the  other  two  methods  although  the  overall  quan¬ 
tification  error  depended  on  the  heterogeneous  arrangement  and 
optical  property  combination  of  heterogeneity  and  background 
calibration  medium  (for  normalization  II).  Furthermore,  it  was 
observed  that  superficially  seated  objects  are  reconstructed 
with  superior  accuracy  than  deep-seated  objects.  Therefore,  the 
observations  in  Figs.  4-7  correspond  to  a  worse  case  scenario. 
This  is  intuitively  correct  since  the  intrinsic  field  measurement 
maps  almost  identically  the  excitation  field  for  the  surface 
fluorochromes  and,  therefore,  the  effects  of  background  het¬ 
erogeneity  cancel  out  when  using  normalization  I.  In  addition, 
the  insensitivity  of  the  normalization  I  and  II  to  unknown  ex¬ 
perimental  and  modeling  inconsistencies,  such  as  the  unknown 
gains  and  attenuations  of  the  system  and  inaccuracies  of  the 
diffusion  approximation  is  a  significant  advantage  of  normal¬ 
ization  methods  leading  to  reduced  experimental  complexity 
and  better  imaging  performance.  The  normalization  III  does  not 
attain  this  advantage  and,  therefore,  reconstructions  become 
sensitive  to  system  calibration  and  accuracy  of  the  theoretical 
model  used.  In  our  case  this  method  resulted  in  compromised 
image  quality  even  when  imaging  the  fluorochrome  tube  alone. 
It  is  possible  that  these  distorted  reconstructions  are  due  to  in¬ 
accurate  modeling  of  the  source  distribution  which  is  assumed 
herein  isotropic,  an  assumption  that  is  not  entirely  consistent 
with  the  exact  photon  propagation  of  the  incident  beam  in  the 
first  few  millimeters  of  propagation.  Since  the  normalization 
I  and  II  are  both  independent  of  the  source  term  they  are  not 
sensitive  to  such  theory-experiment  mismatch. 

Finally,  while  the  current  study  examined  the  effects  of  op¬ 
tical  property  variation,  it  is  expected  that  inhomogeneous  flu¬ 


orescence  distribution  will  further  add  to  the  complexity  of  the 
problem  and  such  investigations  are  currently  considered. 

In  conclusion,  we  have  found  that  inversion  models  based  on 
perturbation  solutions  of  the  diffusion  approximation  can  yield 
robust  imaging  performance  even  at  increased  background  op¬ 
tical  heterogeneity  when  appropriate  data  normalization  is  con¬ 
sidered.  Such  methods  offer  superior  imaging  performance  in 
in  vivo  application  compared  to  nonnormalized  methods  when 
using  linear  inversion  schemes.  This  study  demonstrated  that  the 
use  of  the  normalization  I  (normalized  Bom  ratio)  approach  can 
yield  quantification  accuracy  of  the  order  of  ±10%  as  a  function 
of  background  optical  property  heterogeneity,  especially  when 
calibrations  are  obtained  in  heterogeneous  media.  While  the  ex¬ 
periments  performed  were  modeled  after  the  dimensions,  op¬ 
tical  properties  and  spatially  varying  heterogeneity  seen  in  small 
animal  optical  imaging,  the  normalization  methods  investigated 
in  this  work  could  be  potentially  beneficial  for  clinical  imaging 
as  well,  for  example  breast  optical  tomography. 
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